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Using an approximation method for eigenvalue distribution functions, we study the tempera- 
ture dependence of specific heat of the antiferromagnetic Heisenberg model on the asymmetric 
railroad-trestle lattice. This model contains both the sawtooth-lattice and Majumdar-Ghosh 
models as special cases. Making extrapolations to the thermodynamic limit using finite size 
data up to 28 spins, it is found that specific heat of the Majumdar-Ghosh model has a two- 
peak structure in its temperature dependence and those of systems near the sawtooth-lattice 
point have a three-peak structure. 
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1. Introduction 

Fully frustrated spin systems have been treated ex- 
tensively in both experimental and theoretical investi- 
gations. 1 By the term of fully frustrated spin system, 
we usually mean a system whose classical ground state 
has infinite continuous degeneracies. 2 Some systems have 
been known to have this property. 

A typical example of fully frustrated spin systems is 
the antiferromagnetic Heisenberg model (AFHM) on a 
kagome lattice. Harris et al. studied an ordered ground 
state of the classical kagome AFHM, 3 and pointed out 
that there exist continuous local distortions to the classi- 
cal ground state without changing the total energy. The 
existence of such continuous local distortions leads to the 
ground state degeneracy cx const^ with system size N. 
The quantum dynamics in the spin- 1/2 system makes 
the residual entropy release, and thus the temperature 
dependence of the specific heat shows a two peak struc- 
ture. 4 The unique ground state of the quantum system 
is now believed to be a resonating-valence-bond state. 5 

Spin- 1/2 AFHM on a sawtooth lattice (or A-chain), 
which is one of the simplest one-dimensional fully frus- 
trated spin systems, has been treated by several au- 
thors. Kubo pointed out the twofold-degenerate exact 
dimer-singlet ground state and a double peak structure 
in the specific heat by using the quantum transfer ma- 
trix method. 2 Otsuka applied an approximation method 
for the eigenvalue distribution function (EvDF) to calcu- 
late thermodynamics of AFHM on finite size sawtooth- 
lattices with up to 26 spins under the periodic boundary 
condition. 6 He examined the magnetic-field dependence 
and estimated lower-temperature peak position of the 
specific heat in the thermodynamic limit. The sawtooth 
lattice is expected to be realized in an oxygen doped Cu 
based delafossits YCUO2.5 7 ' 8 an d a molecular quantum 
spin system MO75V20- 9 ' 10 The latter consists of twenty 
spins, and is often called "spin doughnut" because of the 
way of arrangement of S = 1/2 vanadium ions in the 
molecular. 

A lot of theoretical studies have been done for the spin- 



1/2 antiferromagnetic XX Z model on a railroad-trestle 
lattice, because this model shows interesting quantum 
phase transitions and critical phenomenena. 11 Nomura 
and cowarkers determined the quantitative phase di- 
agram by using the level spectroscopy technique. 12-14 
This model contains the Majumdar-Ghosh model as 
a special case. 15 ' 16 The Majumdar-Ghosh model has 
been known to have the twofold-degenerate exact dimer- 
singlet ground state, but it is not a fully frustrated sys- 
tem. 2 

Recently, several authors studied a generalized model, 
where the strength of the leg interactions in the railroad- 
trestle lattice model is extended to be asymmetric. 17-20 
The Hamiltonian is described by 



N 



N 



n = l^ + «Xi {1/2 - (-I)'**} ht i+2 , (i) 



5f 5? + A(Sf Sf + SVS?) (see Fig. 1). The 



with ftAj 

operator Si denotes the quantum spin of size 1/2 on ith 
lattice point, and the periodic boundary condition, i.e., 
Sn+i = Si is imposed. 

The parameter S measures the asymmetry in the leg 
interactions. The ground state phase diagram of this 
model has been investigated by Nakane, Oguchi and one 
of the present authors. 20 Using the level spectroscopy 
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Fig. 1. Schematic representation of the railroad-trestle lattice 
with asymmetric leg interactions. The model reduces to the 
Majumdar-Ghosh (sawtooth-lattice) model at a = 1 and 5 = 
(5 = 1/2). 
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method, they studied how the asymmetry parameter 8 
alters phase boundaries. They found no qualitative effect 
of 8 on the phase boundaries. In particular, the 8 depen- 
dence of the dimer-Neel phase boundary is quite small. 
However, as long as we know, there is no systematic study 
of the thermodynamic properties of this model. 

The main purpose of the present paper is to study the 
thermodynamic properties of the railroad-trestle lattice 
model with asymmetric leg interactions, paying atten- 
tion to the case of a — 1. We set H = H a =i, and con- 
centrate on H below. It is also assumed that < A and 
< 8 < 1/2, and then H has the exact singlet dimer state 
as the ground state. 20 It should be noted that Hs=o is 
the Majumdar-Ghosh model and H$ = i/ 2 is the sawtooth- 
lattice model, as schematically shown in Fig. 1. As men- 
tioned before, the sawtooth-lattice model with 8 = 1/2 
is a fully frustrated system. 2 The thermodynamic prop- 
erties was already studied, and have been known to have 
a two peak structure in the temperature dependence of 
the specific heat. 2,6 The cases of S < 1/2 are not fully 
frustrated, following the previous definition of this term. 
However, the Ising model Hs_a=o has a rather large resid- 
ual entropy even for 8 < 1/2, as described in Sec. 2 in 
detail. This observation motivates us to study the 8 de- 
pendence of the specific heat of the Heisenberg model 

This paper is organized as follows: we study the spe- 
cific heat and entropy of the Ising model Ha=o in § 2, in 
which wc apply the transfer matrix method to this sys- 
tem and the origin of the residual entropy is explained. 
Our calculation process of thermodynamic properties by 
using the EvDF method is explained in § 3. We show 
the calculated results of the Heisenberg model H\ = i in 
§ 4, where we use finite-size data up to 28 sites to make 
an extrapolation to the thermodynamic limit. Finally we 
summarize our results in § 5. 

2. Thermodynamic Properties at the Ising Limit 

In this section, we study the entropy and specific heat 
of the Ising limit of the our model, i?A=o- In order to use 
the transfer matrix method, we rewrite the Hamiltonian 
as the sum of plaquette elements: 

#A=o = 5>(0, (2) 

2 = 1 

where the Ith plaquette Hamiltonian h(l) is denned by 

QZ CZ I QZ QZ 

_ D 2l-l D 2l T 2i+l°2i+2 , Q z qz 



— 8 ) S 2l S 2l+2 , (3) 



+ ( g + 5 J S 2l-l S 2l+l 



p — N/2 is the total number of plaqucttes (or 



and AT, 

triangle elements for the sawtooth-lattice case with 8 
1/2) and assumed a large even number. 

The partition function of the system is given by 



% — IT ^(^Ij-li $2h &2l+li $21+2) 

[S*}1=1 



(4) 



where the sum over {S z } run on all 4 Wp spin configura- 
tions, and 

P(S£i-i> $2h S21+11 $21+2) = e /3h< '^ (5) 

with (3 = 1/T. It is convenient to introduce a 4x4 trans- 
fer matrix P by 

I e -3/3/4 e -/3<5+/4 e /35 + /4 
e /3<5_/4 e /3/4 e /3/4 
e -06-/4 e /3/4 e /3/4 

e /3/4 e /35+/4 e -/3<5+/4 



P = 



\ 



e /3/4 \ 
e -P5-/4 

e -3/3/4 J 



(6) 



(7) 



where 8± = 1 ± 28. Then, we have 

Z = Tr P w p. 

The largest eigenvalue of P is given by 

_ 3e^/ 4 + e- 3 ^/ 4 
^ — z 



+ I v / 5e /3/2 + 4p 5/3 + Ae -8f] + 2e -P/2 + e-3/3/2^) 

Thus we obtain partition function in the limit A^ p — > 00 
as follows: 



An 



Z = X 

which gives the ground state energy per a spin 
1 d 



(9) 



E(T = 0) 



2d(3 



log A 



= -1/8, (10) 



and the residual entropy per a spin 



S(T = 0): 



/3=oo 



ilog3 



for 8 = 1/2 



1 log ii^/I for < 8 < 1/2 



• (11) 



We here comment on the degeneracy factor of the 
ground state manifold which gives the energy in eq. (10). 
For the sawtooth lattice case with 8—1/2, the degener- 
acy factor has been already known. 21 ' 22 It is easy to see 
that the ground state energy E(T = 0) = —1/8 is ob- 
tained if each triangle has two antiferromagnetic bonds 
and one ferromagnetic bond. There are three ways to 
choose the position of a ferromagnetic bond in a trian- 
gle, so that the residual entropy in the first line of eq. (11) 
is obtained. 

Wc turn to the general case of 8 7^ 1/2. It is convenient 
to introduce the terms, "tip spin" and "bottom spin" , as 
shown in Fig. 2. In order to count the degeneracy factor, 
we first choose a bottom spin configuration, 5b- In 5b, 
some bottom spin pairs are in ferromagnetic state and 
the others are in antiferromagnetic state. We assume 5b 
has N-p ferromagnetic bottom spin pairs on ?i, l 2 , • • • , 
ZAr F th triangles. We also define l = In f — N p , taking into 
account the periodic boundary condition. It should be 
noted that the direction of the tip spin on a triangle with 
a ferromagnetic bottom spin pair is fixed to be opposite 
to that of the bottom spins in the ground state manifold 
(see Fig. 2). We now consider how to choose configuration 
of tip spins between the fixed tip spins on the l q and 
Zg_ith triangles, which are shown by the gray circles in 
the right panel in Fig. 2. The tip spin configurations are 
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tip spin 




(a) 



bottom spin 



Fig. 2. Definition of "tip spin" and "bottom spin" (left panel), 
and a schematic explanation of the counting problem of the de- 
generation factor for a fixed bottom spin configurations (right 
panel). The circles represent spins, and the spin direction of the 
open circles is opposite to that of the black circles. The tip spins 
represented by open and black circles are fixed tip spins. The 
exchange coupling surrounded by the gray line determines the 
directions of the gray tip spins between the two fixed tip spins. 



determined so as the tip spin interaction energy between 
the two fixed tip spins to be minimized. The tip spin 
interaction energy contains l q — l q -\ bonds. It is easy to 
observe that the tip spin interaction energy is minimized 
when we put l q — l q -\ — 1 bonds to be antiferromagnetic 
and the remaining one bond to be ferromagnetic. The 
total number of such tip spin configurations is l q — l q -\, 
because of the l q — l q -\ ways to choose one ferromagnetic 
bond. Therefore, for the fixed bottom spin configuration 
5b, the degeneracy factor D(Sb) is given by 



(12) 



Summing up this degeneracy factor over all bottom spin 
configurations, then we obtain 

e NS(T=o) = Y / D(S h ) for 6 < 1/2. (13) 

5 b 

We have numerically checked that this expression repro- 
duces the second line of eq. (11). 

It is instructive to note that D(Sb) is altered to be 
D(S h ) = 2 n p~ Nf for (5 = 1/2, because all configurations 
of iV p — Np tip spins, other than Np fixed tip spins, are 
possible in the ground state manifold. Thus, we obtain 
the first line of eq. (11) again as follows: 



3 JVS(T=0) . 



5 b 



£ ^ 

N F =0 



Cjv„2 



N D -N F 



= 3 



for 6 = 1/2. 



(14) 



We show S(T) and C(T) in Fig. 3 for several val- 
ues of 6. We can confirm in Fig. 3 (a) that S(T) de- 
creases toward to (1/2) In 3 ~ 0.549 for 6 = 1/2 and to 
(1/2) ln[(3 + \/5)/2] ~ 0.481 for the others as decreasing 
temperature. We also find that S(T) for S = 0.4 and 0.3, 
which are close to the sawtooth point S = 0.5, are not 
monotonic: these have a shoulder around T ~ 0.1. This 
behavior of the entropy leads to a low temperature peak 
(or shoulder) in C(T), as seen in Fig. 3(b). 

The present results interest us to study effects of the 
quantum dynamics via spin flip terms. We are going to 
study the Hcisenberg version of this model in the follow- 
ing sections, where it is shown that specific heat of the 
Majumdar- Ghosh model has a two-peak structure in its 
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Fig. 3. Temperature dependence of (a) entropy S(T) and (b) spe- 
cific heat C(T) of the Ising limit for several values of 8. 



temperature dependence and those of systems near the 
sawtooth-lattice point have a three-peak structure. 

3. Numerical Method for an Eigenvalue Distri- 
bution Function and Thermodynamic Prop- 
erties 

In our calculation of specific heat of the Heisenberg 
model, we make an estimation of the eigenvalue distribu- 
tion function (EvDF) with the aid of the Lanczos method 
and a sampling technique, as done by Otsuka. 6 We here 
give a review of this method, and describe our values of 
the parameters used in this method. 

The EvDF p(u) is defined as 



p(uj) = Tr S(uj - H) 

= - -Tr Im.g(u) + ir]), 



(15) 



where H is a Hamiltonian, Q{z) = (z — H)^ 1 , and r\ 
is a real positive infinitesimal, which we approximate as 
2 x 10~ 5 in our numerical calculation. Introducing a basis 
set {(()^\n = 1 <~ N^} in the pth subspace of H, then 
the EvDF can be rewritten in the following form: 



i iV " 

fj. n=l 

The dimension of the subspace, N^, is usually so large 
that we cannot sum over all the basis vectors, which 
leads us to use a sampling technique. In our practical 
calculation, we samples randomly a set of N s vectors, 



Mi 



1 <~ N s }, in such the subspace. Then p^uj) is 
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approximated by 



N S 



7T ' » 



(17) 



71=1 



with the normalization constant c, which is determined 
from the condition J p^(oj)dui = N^. This type of sam- 
pling was used by Imada and Takahashi in their quantum 
transfer Monte Carlo method, 23 where it was discussed 
that enlargement of system size gets N$ to be reduced 
to evaluate physical quantities within a fixed accuracy. 

We turn to the estimation of (^\G(z)\^) = 
Gn n {z). We make a tri-diagonalization of 7i by the Lanc- 

zos method, setting |$^) as the initial vector. Using the 
resultant Lanczos coefficients {ai,f3i\l = 1,2, •••}, the 
matrix element Q% n can be expressed as the following 
continued fraction form: 

= l ~^T2 • (18) 



Z — OL\ — 



OL-2 



102 



It has been known that the convergence of the continued 
fraction form is so fast that we only need to calculate first 
50-200 Lanczos coefficients. 6 Our typical iteration num- 
ber is about 100, for which we checked the convergence 
of calculated quantities using 16-site systems. 

In order to estimate thermodynamic properties, we 
have to calculate the mth moment 



/OO 
u m p(uj)e-P"dLO 
-OO 



(19) 



as a function of j3. It is not efficient to carry out the 
numerical integration over to at each value of (3, because 
p(w) fluctuates much greater than the other factors. Tak- 
ing into account this fact, we first calculate the histogram 

puJi+Su/2 

p(uj)duj, (20) 

-Su/2 

with the energy mesh 8u>. This numerical integration 
should be carried out carefully, because p{u) consists of 
Lorentzians with a very small width of rj. If we choose 
the energy mesh as (35u/2 <C 1, then 

Y 



-P(ui±8u/2) 



-Put 



it 



-@UJi 



(21) 



holds with high accuracy. Therefore, the calculation of 
mth moment reduces to 



(22) 



which give us much efficient way to evaluate the tem- 
perature dependence. Once the temperature dependence 
of the moments is obtained, the specific heat per a spin, 
C(T), is calculated by 



C(T) = 



N 





\Fi(py 


■} 




IMP). 





(23) 



In our calculation, setting our maximum inverse temper- 
ature as /3 max = 100, we adopt fo = 2x 10~ 4 so as to 
be (3Slo/2 < 10~ 2 . 



Finally, we introduce our extrapolation method from 
finite-size data to the thermodynamic limit N — > oo. We 
employ the next equation as the extrapolation function: 



C(T, N) = C(T, N -» oo) 



— (A : Constant) (24) 



where k denotes a positive integer number and later we 
consider about k. 

4. Multi-Peak Structure of the Specific Heat of 
AFHM on the Railroad- Trestle Lattice with 
Asymmetric Leg Interactions 

In this section, we present our numerical results at the 
Heisenberg point A = 1 based on the EvDF method. The 
Hamiltonian treated here is 



N 



Ha=i — J]] Si ■ Si + i 



N/2 

E 

i=l 



£>2i ' £>2(i+l) 



5) So 



'2i+l 



(25) 



Before describing results for < S < 0.5, which is the 
main concern in this paper, we examine our calculations 
at 5 = 0.5, where we can use the existing results. 6 

4-1 Sawtooth- Lattice Point S = 0.5 

In Fig. 4, we show low-temperature specific heat of the 
22-site sawtooth-lattice AFHM together with the unbi- 
ased variance a 2 . Here, N s = 250 is chosen to calculate 
C(T), and we divide the 250 samples into 5 sets of 50 
samples to get the variance, i.e., 

o-\T) = \jy{T)-C i {T)f, (26) 

where Cj (T) represents the specific heat calculated with 
the jth set. We find in Fig. 4 that the sampling error 




Fig. 4. Calculated results of (a) C(T) for the sawtooth-lattice 
AFHM with N = 22, and (b) the unbiased variance <r 2 (T). 
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becomes severe when temperature decreases and C(T) 
varies rapidly. As described later, sampling errors at tem- 
peratures below the peak temperature, T p ~ 0.06, poses 
an obstacle to our extrapolation procedure at low tem- 
peratures. 

We now study the dependence of C{T) on the sys- 
tem size TV (even number). We show C(T) in Fig. 5, 
classifying by TV/2 being an even or odd number to ob- 
tain smooth evolution on the system size. If we don't 
make such a classification, for example, the peak height 
as a function of TV shows a nonmonotonic behavior. In 
order to understand this behavior, we focus on the bot- 
tom structure of TV-site sawtooth lattice. A TV-site saw- 
tooth lattice has a bottom shaped polygon which has 
TV/2 angles. An example is given in Fig. 6, where we 
show a sawtooth lattice with N = 10 whose bottom 
structure is a pentagon. Let us assume that the cou- 
pling for the bonds represented by the dashed line in 
Fig. 6 can be omitted. Then, only the bottom polygon 
contributes to specific heat. The ground-state of the pen- 
tagon in Fig. 6(b) consists of one ferromagnetic bond and 
four antiferromagnetic bonds. We call this situation "1/5 
frustration" . On the other hand, if the bottom structure 
is square, hexagon and so on, then there exists no frustra- 
tion. Thus, if N/2 is odd number, the bottom polygon in 



(a) 



(b) 




Fig. 5. Finite-size and extrapolated results of C(T) for the 
sawtooth-lattice AFHM, where N/2 is even in (a) and odd in 
(b). The offset in the horizontal axis by 1/N is applied for the 
A^-spin results. 



© ,© © ® © 
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(a) 



Fig. 6. Graphical representations of (a) 10-site sawtooth lattice 
and (b) its bottom structure. 




the TV-site sawtooth system has 2/TV frustration. In this 
argument, we have omitted the dashed-line bonds. How- 
ever, the nonmonotonic behavior on TV in our calculated 
data seems to indicate that such a simple picture holds 
even when the omitted couplings are again taken into 
account. In our extrapolation procedure, the two data 
sets are treated independently and the consistency be- 
tween the two extrapolated results is used to check our 
procedure. 

We turn to the index k defined in (24). To estimate 
an adequate value of k, we examine the specific heat as 
a function of TV for fixed temperatures, T = 0.15 and 
0.125, at which sampling error is not so large that the 
size dependence is not smeared. In Fig. 7 the specific 
heat Ct(TV) against 1/TV 3 is shown. We find a linear 
behavior for the large systems, and thus we employ the 
extrapolation function (24) with k = 3 in the case of 
sawtooth lattice. The extrapolation results for TV — > oo 
are shown in Fig. 5 as the open circles. The extrapolation 
by the even-TV/2 data set in Fig. 5(a) is consistent with 
that by the odd- TV/2 one, as expected. 



S. 0.18 



A 7=0.15 
V 7=0.125 



4 

1W 3 



Fig. 7. Plot of C T (N) vs 1/N 3 at T = 0.15 and 0.125, where 
<5 = 0.5. 



Otsuka calculated the specific heat by the same 
method for odd- TV/2 systems up to TV = 26. 6 Our cal- 
culated results are in good agreement with his results at 
temperatures higher than the peak temperature T p . Al- 
though both results are also consistent with each other 
within the sampling error at T < T p , there exists small 
deviation in average values. Magnitude of the deviation 
is order of the marker size in Fig. 5. This small deviation 
gets our extrapolation result to be somewhat larger than 
Otsuka's extrapolation at low temperatures. Expanding 
the sampling number may resolve this problem. How- 
ever, we don't pursue this problem, because it is time- 
consuming and docs not change our qualitative conclu- 
sions. 

4-2 Majumdar- Ghosh Point 5 = 

We here study the Majumdar- Ghosh model with 5 = 
0. This model is well known for the exact singlet-dimer 
ground state. 15 ' 16 However, as long as the present au- 
thors know, the thermodynamic property has never been 
examined. 
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We show the temperature dependence of the low tem- 
perature specific heat for finite size systems with TV = 16, 
20, 24 and 28 in Fig. 8. (For the N = 16 system, the ex- 
act result is also shown as the open circles, which are in 
good agreement with the EvDF result.) We find in Fig. 8 
that a shoulder appears around T = 0.1 as the system 
size is enlarged. This size dependence make us antici- 
pate the appearance of a low-temperature peak in the 
thermodynamic limit. In fact, following the extrapola- 
tion procedure explained in the previous section, we ob- 
tain a low-temperature peak in our extrapolation shown 
by the solid line in Fig. 8. The peak temperature is esti- 
mated as Tp =0 w 0.09, which is higher than that of the 

sawtooth-lattice model, Tp -1 ^ 2 w 0.06. 

In Fig. 9, we show the specific heat of N — ► oo for 
a wider temperature range. It is found that a high- 
temperature peak exists at T w 0.3, which is almost 




T 



Fig. 8. Low-tcmpcraturc specific heat of the Majumdar-Ghosh 
model with S = 0.0 for finite-size systems and the thermody- 
namic limit. The open circles represent the exact specific heat 
for the 16-site system. 



same as the peak temperature of the Ising limit (see 
Fig. 3). This result indicates that the spin flip terms in 
the Majumdar-Ghosh model get the residual entropy of 
Ising limit to release and lead to the low-temperature 
peak at T w 0.09. 

The two-peak structure of the specific heat for the 
Majumdar-Ghosh model is one of the main conclusions 
in this paper. In order to check this further, we carry 
out a preliminary calculation on the 30-site Majumdar- 
Ghosh model. The result is shown in Fig. 9 as the dashed 
line, together with sampling error. Although it has large 
sampling error due to smallness of the sampling number, 
it suggests that the 30-site Majumdar-Ghosh model has 
a two-peak structure in the temperature dependence of 
the specific heat, which gives a support to our conclusion. 

4-3 Evolution of the Specific Heat as a Function of 8 

Finally, we study how C(T) depends on the asymmetry 
parameter 5. In Fig. 10, we show C(T) of the 24- and 26- 
site systems for S = 0.0, 0.1, • • • , 0.5. In both systems, 
C(T) for S = 0.0, 0.1 and 0.2 has a high-temperature 
peak at T ~ 0.3 and a shoulder around T = 0.1. For 
S = 0.3, the high-temperature peak and the shoulder 
merge into a terracing. For 5 = 0.4 and 0.5, C(T) shows a 
two-peak structure clearly. In addition, a small shoulder 
at T ~ 0.1 appears. 

Making extrapolation to the thermodynamic limit, we 
obtain the results shown in Fig. 11. Comparing with the 
finite size data, we find the low-temperature shoulder 
for S = 0.0, 0.1 and 0.2 changes to a low-temperature 
peak. Also, the low-temperature side edge of the terrac- 
ing for S = 0.3 is separated from it and turns into a 
low-temperature peak. For S = 0.4, the small shoulder 
at T ~ 0.1 changes to a small peak, and thus C(T) for 
5 = 0.4 has three peaks. The small shoulder for 5 = 0.5 
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Fig. 9. (a) Specific heat for the whole temperature range for the 

Majumdar-Ghosh model with N — > oo and our preliminary result 
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Fig. 11. Estimations of C(T) in the thermodynamic limit for the 
AFHM on the asymmetric railroad-trestle lattice. 



becomes clearer, but it does not turn to a peak. Thus, we 
conclude that the present model has a two-peak structure 
in the specific heat, and an additional mid-temperature 
peak appears for S ~ 0.4. 

To understand the overall behavior of the specific heat, 
it is helpful to make a comparison with the result of the 
Ising limit. As for <5-dependence of the high-temperature 
peak, the Heisenberg system has similar features to the 
Ising limit as seen in Figs. 3 and 11: for the both cases, 
the high-temperature peak position as a function of S is 
almost same, and the peak height is decreased as 5 in- 
creasing. The lowest-temperature peak in Fig. 11, which 
originates in the spin flip terms, grows when S increases. 
This observation can be interpreted in terms of the Ising 
limit: we find in Fig. 3(a) that, at low temperatures, the 
entropy is larger for larger 5, and thus the spin flip term 
can potentially lead to larger lowest-temperature peak 
for larger i5. Also, in terms of the Ising limit, the three- 
peak structure for S = 0.4 in Fig. 11 is plausible, because 
the specific heat has already two peaks in the Ising limit. 
Finally, we notice in Fig. 11 that C(T) becomes to be 
more sensitive to S as S increasing. We can find the same 
tendency in Fig. 3 for the Ising limit. 

5. Summary 

In this paper, we consider thermodynamic proper- 
ties of the Ising and Heisenberg antifcrromagncts on a 
railroad-trestle lattice with asymmetric leg interactions. 
At first we have investigated entropy structure of the 
sawtooth-lattice and railroad-trestle-lattice Ising antifcr- 
romagnets, where we have found that the railroad-trestle- 
lattice model has less residual entropy than the sawtooth- 
lattice model because of difference in total number of 
possible tip-spin configurations in the ground state man- 
ifold. Next we have calculated thermodynamics of the 
Heisenberg antiferromagnets by using an eigenvalue dis- 
tribution function. Our obtained results are summarized 
as follows. (1) We have found a finite-size effect: we can 
classify finite-size results of C(T) according to TV/2 be- 
ing even or odd. This tendency becomes to be more 
pronounced for systems with larger asymmetry. (2) In 



the thermodynamic limit, our calculated results indi- 
cate that the Heisenberg antiferromagnet on the asym- 
metric railroad-trestle lattice has two or more peaks in 
the temperature dependence of the specific heat. Espe- 
cially, the specific heat has a three peak structure for 
8 ~ 0.4, and a two peak structure even at 5 = 0. For 
(5 = 0, the Majumdar-Ghosh model, we estimate the 
lower-temperature peak position as T p ss 0.09. We have 
also pointed out that the 5 dependence of specific heat 
of the Heisenberg model is possible to be interpreted in 
terms of the Ising limit. 

As for finite size data themselves, we have not ob- 
served multi-peak structure in C(T) at S < 0.3 within 
the present calculations for TV < 28, though the size de- 
pendence of C(T) exhibits a sign of the multi-peak struc- 
ture. As a future problem, it is interesting to study the 
specific heat of larger systems than N = 28 by using the 
density matrix rcnormalization group method, 24 which 
provides a critical check of the conclusion obtained in 
this paper. 
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